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1 Introduction 

It is well established to view the integer quantum Hall effect (QHE) as a 
sequence of quantum phase transitions associated with critical points that 
separate energy regions of localised states where the Hall-conductivity a xy 
is quantised in integer units of e 2 /h (see, e.g., |1I2|). Simultaneously, the 
longitudinal conductivity a xx gets unmeasurable small in the limit of van- 
ishing temperature and zero frequency. To check the inherent consequences 
of this theoretical picture, various experiments have been devised to investi- 
gate those properties that should occur near the critical energies E n assigned 
to the critical points. For example, due to the divergence of the localisation 
length oc \E — E n \~>*, the width A of the longitudinal conductivity 

peaks emerging at the transitions is expected to exhibit power-law scaling 
with respect to temperature, system size, or an externally applied frequency. 

High frequency Hall-conductivity experiments, initially aimed at resolving 
the problem of the so-called low frequency breakdown of the QHE apparently 
observed at ~ 1 MHz, were successfully carried out at microwave frequencies 
(~ 33 MHz) The longitudinal ac conductivity was also studied to obtain 
some information about localisation and the formation of Hall plateaus in the 
frequency ranges 100 Hz to 20 kHz 0] and 50-600 MHz 0. Later, frequency 
dependent transport has been investigated also in the Gigahertz frequency 
range below 15 GHz ©El and above 30 GHz |8I0I1()| . 

Dynamical scaling has been studied in several experiments, some of which 
show indeed power law scaling of the <j xx (uj) peak width as expected, A ~ 
w K |filllll2j . whereas others do not [Jj. The exponent k = (^z)^ 1 contains 
both the critical exponent fj, of the localisation length and the dynamical 
exponent z which relates energy and length scales, E ~ L~ z . The value 
of /i = 2.35 ± 0.03 is well known from numerical calculations |13I14| . and 
it also coincides with the outcome of a finite size scaling experiment |15| . 
However, it is presently only accepted as true that z = 2 for non-interacting 
particles, and z = 1 if Coulomb electron-electron interactions are present 
|16I17I18I1HI20| . Therefore, a theoretical description of the ac conductivities 
would clearly contribute to a better understanding of dynamical scaling at 
quantum critical points. 
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Legal metrology represents a second area where a better knowledge of 
frequency dependent transport is highly desirable because the ac quantum 
Hall effect is applied for the realization and dissemination of the impedance 
standard and the unit of capacitance, the farad. At the moment the achieved 
relative uncertainty at a frequency of 1 kHz is of the order 10 -7 which still is 
at least one order of magnitude to large |21I22I23I24| . It is unclear whether the 
observed deviations from the quantised dc value are due to external influences 
like capacitive and inductive couplings caused by the leads and contacts. Al- 
ternatively, the measured frequency effects that make exact quantisation im- 
possible could be already inherent in an ideal non-interacting two-dimensional 
electron gas in the presence of disorder and a perpendicular magnetic field. Of 
course, applying a finite frequency u> will lead to a finite a xx (u) and in turn 
will influence the quantisation of a xv (u>), but it remains to be investigated 
how large the deviation will be. 

Theoretical studies of the ac conductivity in quantum Hall systems started 
in 1985 [2SJ when it was shown within a semiclassical percolation theory that 
for finite frequencies the longitudinal conductivity is not zero, thus influencing 
the quantisation of the Hall-conductivity. The quantum mechanical problem 
of non-interacting electrons in a 2d disordered system in the presence of 
a strong perpendicular magnetic field B was tackled by Apel (2E] using a 
variational method. Applying an instanton approximation and confining to 
the high field limit, i.e., restricting to the lowest Landau level (LLL), an 
analytical solution for the real part of the frequency dependent longitudinal 
conductivity could be presented, a xx (u>) ex w 2 ln(l/w 2 ), a result that should 
hold if the Fermi energy lies deep down in the lower tail of the LLL. 

Generalising the above result, both real and imaginary parts of the fre- 
quency dependent conductivities were obtained in a sequence of papers by 
Viehweger and Efetov 27 28 29 . The Kubo conductivities were determined 
by calculating the functional integrals in super-symmetric representation near 
non-trivial saddle points. Still, for the final results the Fermi energy was re- 
stricted to lie within the energy range of localised states in the lowest tail of 
the lowest Landau band and, therefore, no proposition for the critical regions 
at half fillings could be given. The longitudinal conductivity was found to be 

a xx {uj) = cuj 2 ln b (l/w 2 ) - ie 2 Lu2l 2 B g{E F ) , (1) 

with the density of states g(Ep), and two unspecified constants b and c |2*7) . 
The real part of the Hall-conductivity in the same limit was proposed as 

a yx {u) = e 2 /h (huj/T) 2 8Trl 2 B n e , (2) 

where Ib = WK/JeB) is the magnetic length, r 2 — X/2nl 2 B is the second mo- 
ment of the white noise disorder potential distribution with disorder strength 
A, and n e denotes the electron density. The deviation due to frequency of the 
Hall-conductivity from its quantised dc plateau value can be perceived from 
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an approximate expression proposed by Viehweger and Efetov |29| for the 
2nd plateau, e.g., for filling factor i/f»2, 

= e 2 /h [ x _ ^ - (^ff (2 - „)] . (3) 

Again, r is a measure of the disorder strength describing the width of the 
disorder broadened Landau band, and uj c = eB/m e is the cyclotron frequency 
with electron mass m e . According to ©, due to frequency a deviation from 
the quantised value becomes apparent even for integer filling. 

Before reviewing the attempts which applied numerical methods to over- 
come the limitations that had to be conceded in connexion with the position 
of E-p in the analytical work, and to check the permissiveness of the approx- 
imations made, it is appropriate here to mention a result for the hopping 
regime. Polyakov and Shklovskii [3U| obtained for the dissipative part of the 
ac conductivity a relation which, in contrast to Q, is linear in frequency oj, 

fta xx {u;) = Ke& . (4) 

This expression has recently been used to successfully describe experimental 
data |31|, Here, £ is the localisation length, e the dielectric constant, and the 
pre-factor is K = 1/6 in the limit hu ^> k^,T. Also, the frequency scaling of 
the peak width A ~ u> K was proposed within the same hopping model |30| . 

Turning now to the numerical approaches which were started by Gammcl 
and Brenig who considered the low frequency anomalies and the finite size 
scaling of the real part of the conductivity peak a^ x (uj,Ly) at the critical 
point of the lowest Landau band |S2I ■ For these purposes the authors utilised 
the random Landau model in the high field limit (lowest Landau band only) 
|1I13| and generalised MacKinnon's recursive Green function method for 
the evaluation of the real part of the dynamical conductivity. In contrast to 
the conventional quadratic Drude-like behaviour the peak value decayed lin- 
early with frequency, cr^. x (ui) = a xx — const. \u>\, which was attributed to the 
long time tails in the velocity correlations which were observed also in a semi- 
classical model [34135) . The range of this unusual linear frequency dependence 
varied with the spatial correlation length of the disorder potentials. A second 
result concerns the scaling of a^ x (u> , L y ) at low frequencies as a function of the 
system width L y , cr^fwL 2 ) oc (wi 2 ) 2 "''/ 2 , where T) = 2-D(2) = 0.36 ±0.06 
|3*2~] is related to the multi-fractal wave functions |36l37l38l3l?| . and to the 
anomalous diffusion at the critical point with rj ~ 0.38 0D] an d D{2) ~ 1.62 

im 

The frequency scaling of the a xx (u>) peak width was considered numeri- 
cally for the first time in a paper by Avishai and Luck |42j . Using a contin- 
uum model with spatially correlated Gaussian disorder potentials placed on 
a square lattice, which then was diagonalised within the subspace of func- 
tions pertaining to the lowest Landau level, the real part of the dissipative 
conductivity was evaluated from the Kubo formula involving matrix elements 
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of the velocity of the guiding centres This is always necessary in single 
band approximations because of the vanishing of the current matrix elements 
between states belonging to the same Landau level. As a result, a broadening 
of the conductivity peak was observed and from a finite size scaling analysis 
a dynamical exponent z = 1.19 ± 0.13 could be extracted using v = 2.33 
from ^3]- This is rather startling because it is firmly believed that for non- 
interacting systems z equals the Euclidean dimension of space which gives 
z = d = 2 in the QHE case. 

A different theoretical approach for the low frequency behaviour of the 
&xx(w) peak has been pursued by Jug and Ziegler |44j who studied a Dirac 
fermion model with an inhomogeneous mass |45| applying a non-perturbative 
calculation. This model leads to a non-zero density of states and to a finite 
bandwidth of extended states near the centre of the Landau band . There- 
fore, the ac conductivity as well as its peak width do neither show power-law 
behaviour nor do they vanish in the limit lo — > 0. This latter feature of the 
model has been asserted to explain the linear frequency dependence and the 
finite intercept at u> = observed experimentally for the width of the conduc- 
tivity peak by Balaban et al. |7|, but, up to now, there is no other experiment 
showing such a peculiar behaviour. 

2 Preliminary considerations — basic relations 

In the usual experiments on two-dimensional systems a current I x {lo) is driven 
through the sample of length L x and width L y . The voltage drop along the 
current direction, U x (co), and that across the sample, U v (uj), are measured 
from which the Hall-resistance Rh(uj) = U y {uj) / I x {tu) = p xy (uj) and the 
longitudinal resistance R x (to) — U x (us) / I x {uj) = p xx (iS)L x j L y are obtained, 
where p xy and p xx denote the respective resistivities. To compare with the 
theoretically calculated conductivities one has to use the relations in the 
following, only in Corbino samples <j xx (uj) can be experimentally detected 
directly. 

The total current through a cross-section, I x (uj) — J^ v j x (uj,r)dy, is de- 
termined by the local current density j x (tv, r) which constitutes the response 
to the applied electric field 



The nonlocal conductivity tensor is particularly important in phase-coherent 
mesoscopic samples. Usually, for the investigation of the measured macro- 
scopic conductivity tensor one is not interested in its spatial dependence. 
Therefore, one relies on a local approximation and considers the electric field 
to be effectively constant. This leads to Ohm's law j = aE from which the 
resistance components are simply given by inverting the conductivity tensor 




u&{x,y} 



(5) 
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Pxx Pxy \ 1 | &yy &yx \ j-g-j 

^Pyx Pyy ) <J xx (Jyy — O xy <J yx \ a xy & xx ) 

For an isotropic system we have @xx — & yy and u yx = —<J X y which in case of 
zero frequency gives the well known relations 

_ ®xx _ O xy 

Pxx = — 5 ; n~ 7 Pxy — — 5 j 5— ■ (,') 

(7 + (7 (7 + (7 

^za; 1 xy w xx 1 

From experiment one knows that whenever p xy is quantised p xx gets immea- 
surable small which in turn means that a xx — ► and p xy = l/cr yx . Therefore, 
to make this happen one normally concludes that the corresponding electronic 
states have to be localised. 

In the presence of frequency this argument no longer holds because elec- 
trons in localised states do respond to an applied time dependent electric field 
giving rise to an alternating current. Also, both real and imaginary parts have 
to be considered now 

a xx (u) = a xx {u) + \ct i xx {u), a xy (u) = a xy {uj) + ia^M . (8) 

Assuming an isotropic system, the respective tensor components of the ac 
resistivity p uv (lo) — p^ v {lo) + ip^ v (ui) with u,v G {x,y} can be written as 

R , )= Vyui^lx + fojy) + ^Li^xx °lx + °xy ^xy) (q) 

Puv[U) {5al x + 6a 2 xy ) 2 + 4(a« a xx + a% a xy ) 2 1 > 

PuvW (8a 2 xx +Sa 2 y ) 2 +A(a^ x a xx + a^al y ) 2 [ > 



with the abbreviations 8a 2 xx = (crfj 2 - (cr^) 2 and 5a 2 xy = (a^ y ) 2 - (a 1 ^) 



2 



3 Model and Transport theory 

We describe the dynamics of non-interacting particles moving within a two- 
dimensional plane in the presence of a perpendicular magnetic held and ran- 
dom electrostatic disorder potentials by a lattice model with Hamiltonian 

ff = 5> r |r><r|- £ V rr ,\r)(r'\. (11) 

r < r ^r'> 

The random disorder potentials associated with the lattice sites are denoted 
by w r with probability density distribution P(w r ) — 1/W within the interval 
[—W/2,W/2], where W is the disorder strength, and the \r) are the lattice 
base vectors. The transfer terms 

r' 

V rr . =Vcxp(-ie 2 /h J A(I)dZ), (12) 
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Fig. 1. Disorder averaged imaginary part of a% x {E,u,£) at energy E/V — —3.35, 
frequency hio/V = 0.001, and e/V = 0.0008 as a function of sample length L x = 
Na. The system width is L y /a — 32 and the number of realisations amounts to 29 
for Lx/a < 5 ■ 10 6 and to 8 for larger lengths 

which connect only nearest neighbours on the lattice, contain the influence 
of the applied magnetic field via the vector potential A(r) = (0, Bx, 0) in 
their phase factors. V and the lattice constant a define the units of energy 
and length, respectively. 

The electrical transport is calculated within linear response theory us- 
ing the Kubo formula which allows to determine the time dependent linear 
conductivity from the current matrix elements of the unperturbed system 



where f(E) = (exp[(_E — E-F)/(k^,T)] — 1) 1 is the Fermi function. The area 
of the system is fl = L x L yi v u — i/Ti[H,u] signifies the u-component of the 
velocity operator, and S(E + huj-H)= i/ {2ir)(G' JJ ' + {E) - G U -{E)), where 
G u,± (E) = ((E + hui ± lt)I — H)^ 1 is the resolvent with imaginary frequency 
ie and unit matrix I. 

For finite systems at temperature T — OK, ensuring the correct order of 
limits for size fl and imaginary frequency ie, one gets with 7 = lo + 2ie /Ti 



a uv (E F ,T,uj) = — 




(13) 
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Fig. 2. The variance of {Q a* x (E , u , e)} versus system length L x calculated for the 
averages shown in Fig. exhibiting an empirical power-law ~ L x ' 5 



lim lim 



£—►0 Q— too h QTiui 

pEf — Klo 



pEp 

(E,uj,e,L x ,Ly)dE 

J Ep—huj 



/n, F—nu) 
-oc 



(15) 



One can show that the second integral with the limits (-co, Ep — tko) does 
not contribute to the real part of a uv (E-p,u>) because the kernel is identical 
zero, but we were not able to proof the same also for the imaginary part. 
Therefore, using the recursive Green function method explained in the next 
section, we numerically studied 3cr*^(i?, lo, e, L x , L y ) and found it to become 
very small only after disorder averaging. As an example we show in Fig. ^ 
the dependence of {Qa* x ), averaged over up to 29 realisations, on the length 
of the system for a particular energy E/V — —3.35 and width L y /a = 32. 
Also the variance of (a* x (E, u>, e)) gets smaller with system length follow- 
ing an empirical power law ~ (L x /a)~ ' 5 (see Fig. |2J. In what follows we 
neglect the second integral for the calculation of the imaginary parts of the 
conductivities and assume that only the contribution of the first one with 
limits (E-p — Hu),Ef) matters. Of course, on has to be particularly careful 
even if the kernel is very small because with increasing disorder strength the 
energy range that contributes to the integral (finite density of states) tends 
to infinity. Therefore, a rigorous proof for the vanishing of this integral kernel 
is highly desirable. 



4 Recursive Green function method 



A very efficient method for the numerical investigation of large disordered 
chains, strips and bars that are assembled by successively adding on slice at a 
time has been pioneered by MacKinnon |47j . This iterative technique relies on 
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Fig. 3. The real and imaginary parts of a xx (E, u>) as a function of energy and 
frequency huj/V = 2-l(T 4 (■), 5-l(T 4 (▼), MO -3 (♦), 2-10" 3 (a). For comparison, 
the cyclotron frequency is hu> c ~ 1.57 V for ub = 1/8 which was chosen for the 
magnetic flux density. The disorder strength is W/V = 1 and the maximal system 
width amounts to L y /a — 96 with periodic boundary conditions applied 

the property that the Hamiltonian H^ N+1 ^ of a lattice system containing iV+1 
slices, each a lattice constant a apart, can be decomposed into parts that de- 
scribe the system containing N slices, H^ N \ the next slice added, iTjv+i jv+ij 
and a term that connects the last slice to the rest, H' N = if^iv+i + fljv+i.iV' 
Then the corresponding resolvent is formally equivalent to the Dyson equa- 
tion G = Gq + GoVG where the 'unperturbed' Go represents the direct sum 
of ifW and Hn+i,n+i, an d V corresponds to the 'interaction' H' N . The es- 
sential advantage of this method is the fact that, for a fixed width L y , the 
system size is increased in length adding slice by slice, whereas the size of the 
matrices to be dealt with numerically remains the same (47148) . 

A number of physical quantities like localisation length (4915015 ll5"2"] . den- 
sity of states (51153) . and some dc transport coefficients (33154155] have been 
calculated by this technique over the years. Also, this method was imple- 
mented for the evaluation of the real part of the ac conductivity in Id |56I57| 
and 2d systems 32 58 59 . Further efforts to included also the real and imagi- 
nary parts of the Hall- and the imaginary part of the longitudinal conductivity 
in quantum Hall systems were also successfully accomplished (61)161162] , 

The iteration equations of the resolvent matrix acting on the subspace of 
such slices with indices i, j < N in the N-th iteration step can be written as 

u i,j — u i,j ' u i,N tlN,N+l <J-jV-fi,Ar + i U N+1,N ^N,j 
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Fig. 4. The imaginary part of a xx (E,u}) as a function of the real part. Data are 
taken from Fig. [3] 
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The calculation of the ac conductivities starts with the Kubo formula H15JI 
by setting up a recursion equation for fixed energy E, width L y = Ma, and 
imaginary frequency e, which, e.g., for the longitudinal component reads 



G 



Hn+i,nG^^ NS> Hn,n+i] 1 



(16) 



a xx (E,w,e,N) 
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(huj) 2 x l Gt- + x ] G+ + 2isS ij x 2 i (G[ 



0J,+ 

ii 



G 



ji 1 



(17) 



The iteration equation for adding a new slice is given by 



with R 



E atR A t i i DxjR 



'N 1L N+l-' 
i4 



DnRn + iD%Rn+i 



(18) 



JV+l.JV+l ' 



and a set of auxiliary quantities as defined in the 
appendix. The coupled iteration equations and the auxiliary quantities are 
evaluated numerically, the starting values are set to be zero. In addition, 
coordinate translations are required in each iteration step to keep the origin 
xn + i = which then guarantees the numerical stability. 



5 Longitudinal conductivity cr xx (E,u)) 

In this section we present our numerical results of the longitudinal conduc- 
tivity as a function of frequency for various positions of the Fermi energy 
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huj/v 

Fig. 5. The real (•) and imaginary (■) part of a xx (Ep,u) in units of e 2 /h at 
E-p/V = —3.5 as a function of frequency. Further parameters are the disorder 
strength W/V = 1.0, system width L y /a = 32 and e/V = 0.0004 

"within the lowest Landau band. The real and imaginary parts of a xx (E-p,Lu) 
were calculated for several frequencies, but, for the sake of legibility, only 
four of them are plotted in Fig. [3] versus energy. While the real part exhibits 
a positive single Gaussian-like peak with maximum « 0.5 e 2 /h at the critical 
energy, the imaginary part, which is negative almost everywhere, has a dou- 
ble structure and vanishes near the critical point. 3 cr xx (Ep, lu) almost looks 
like the modulus of the derivative of the real part with respect to energy. 
Plotting 3 <j X x(Ef, lu) as a function of $la xx (EF,Lu) (see Fig.0J we obtain for 
frequencies lu < lu* a single, approximately semi-circular curve that, up to 
a minus sign, closely resembles the experimental results of Hohls et al. \M\ . 
However, for larger lu our data points deviate from a single curve. 

5.1 Frequency dependence of real and imaginary parts 

The behaviour of the real and imaginary part of the longitudinal ac conduc- 
tivity in the lower tail of the lowest Landau band (E/V — —3.5) is shown 
in Fig. |SJ We find for the imaginary part a linear frequency dependence for 
small lu which is, apart from a minus sign, in accordance with [57]. The 
real part can nicely be fitted to cx lo 2 log 2 [V / \Hlu)] 2 in conformance with Q 
and b — 2, but disagrees with the findings in where b = 1 was proposed. 

5.2 Behaviour of the maximum of er xx (u)) 

The frequency dependence of the real part of the longitudinal conductivity 
peak value was already investigated in [22| where for long-range correlated 
disorder potentials a non-Drude-like decay was observed. We obtained a sim- 
ilar behaviour also for spatially uncorrelated disorder potentials in a lattice 
model jSJ. In Fig. the difference |5R <t xx (lo) — ^^(O)! is plotted versus fre- 
quency in a double logarithmic plot from which a linear relation can be 
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Fig. 6. The non-Drude decrease of the peak value of the longitudinal conductivity 
as a function of frequency. A linear behaviour of |St (T xx (uj) — <?xx(0)\ is clearly 
observed with a% x (Q) / (e 2 / h) = 0.512 using the following parameters E/V = —3.29, 
W/V = 0.1, L y /a = 32, e/V = 0.0004 



discerned. A linear increase with frequency was found for the imaginary part 
of the longitudinal conductivity at the critical point as well [HI]. The stan- 
dard explanation for the non-Drude behaviour in terms of long time tails in 
the velocity correlations, which were shown to exist in a QHE system |32I35| . 
seems not to be adequate in our case. For the uncorrelated disorder potentials 
considered here, it is not clear whether the picture of electron motion along 
equipotential lines, a basic ingredient for the arguments in |32j . is appropri- 
ate. 



5.3 Scaling of the cr xx (u)) peak width 

The scaling of the width of the conductivity peaks with frequency is shown 
in Fig. \7\ where both the <j xx {uj) peak width expressed in energy and, due to 
the knowledge of the density of states, in filling factor are shown to follow 
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Fig. 7. Frequency scaling of the a xx {uj) peak width. The width in energy AE (□) 
and the width in filling factor Av (O) show a power-law A ~ uj k with an exponent 
k = (fJ-z)' 1 = 0.21 
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Fig. 8. Frequency dependent Hall-conductivity: The figure on the l.h.s shows a xy (oj) 
versus filling factor for two system sizes L y /a — 64 (o) and L y /a — 128 (□). 
On the r.h.s., the relative deviation of the Hall-conductivity from the dc value, 
a xy {oj)/cr X y(0) — 1, is plotted as a function of frequency for v = 0.88 



a power-law A ~ uj k with n — (/xz) _1 = 0.21 jHOj. Taking fi = 2.35 from 
[Xl) we get z = 2.026 close to what is expected for non-interacting electrons. 
Therefore, the result reported in 03| seems to be doubtful. However, spatial 
correlations in the disorder potentials as considered in 142} may influence the 
outcome. Alternatively, one could fix z — 2 and obtain /x = 2.38 in close 
agreement with the results from numerical calculations of the scaling of the 
static conductivity or the localisation length |14|. 

The experimentally observed values k ~ 0.42 [H] and k = 0.5 ± 0.1 ^1 
are larger by a factor of about 2. This is usually attributed to the influence of 
electron-electron interactions (z = 1) which were neglected in the numerical 
investigations. 



6 Frequency dependent Hall-conductivity 

The Hall-conductivity due to an external time dependent electric field as a 
function of filling factor v is shown in Fig. [HI for system widths L y /a = 64 
and Ly/a = 128, respectively. While a xy {u)) has already converged for E F 
lying in the upper tail of the lowest Landau band a pronounced shift can be 
seen in the lower tail of the next Landau band. For v = 1.3 a system width 
of at least L y /a = 192 was necessary for <t xv (uj) to converge. This behaviour 
originates in the exponential increase of the localisation length with increas- 
ing Landau band index. Due to the applied frequency Huj/V — 0.008 the a xy 
plateau is not flat, but rather has a parabola shape near the minimum at 
v = 1.0, similar to what has been observed in experiment 24 . An example 
of the deviation of <r xy (uj) from its quantised dc value is shown on the right 
hand side of Fig. [S] where a xy (uj)/a xy (0) — 1 is plotted versus frequency for 
filling factor v = 0.88. A power-law curve ~ w 5 can be fitted to the data 
points. Using this empirical relation, we find a relative deviation of the order 
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of 5 • 10 -6 when extrapolated down to to 1 kHz, the frequency usually applied 
in metrological experiments 21 22 23 24 . Therefore, there is no quantisation 
in the neighbourhood of integer filling even in an ideal 2d electron gas with- 
out contacts, external leads, and other experimental imperfections. Recent 
calculations, however, show that this deviation can be considerably reduced, 
even below 1 • 10 -8 , if spatially correlated disorder potentials are considered 
in the model |64j . 

7 Conclusions 

The frequency dependent electrical transport in integer quantum Hall sys- 
tems has been reviewed and the various theoretical developments have been 
presented. Starting from a linear response expression a method has been 
demonstrated which is well suited for the numerical evaluation of the real 
and imaginary parts of both the time dependent longitudinal and the Hall- 
conductivity. In contrast to the analytical approaches, no further approxi- 
mations or restrictions such as the position of the Fermi energy have to be 
considered. 

We discussed recent numerical results in some detail with particular em- 
phasis placed on the frequency scaling of the peak width of the longitudinal 
conductivity emerging at the quantum critical points, and on the quantisa- 
tion of the ac Hall-conductivity at the plateau. As expected, the latter was 
found to depend on the applied frequency. The extrapolation of our calcu- 
lations down to low frequencies resulted in a relative deviation of 5 ■ 10 -6 
at ~ 1 kHz when spatially uncorrelated disorder potentials are considered. 
Disorder potentials with spatial correlations, likely to exist in real samples, 
will probably reduce this pronounced frequency effect. 

Our result for the frequency dependence of the <r xx peak width showed 
power-law scaling, A ~ u K , where k = (fJ-z)^ 1 = 0.21 as expected for non- 
interacting electrons. Therefore, electron-electron interactions have presum- 
ably to be taken into account to explain the experimentally observed k ks 0.5 
Also, the influence of the spatial correlation of the disorder potentials may 
influence the value of k. 

The frequency dependences of the real and imaginary parts of the longitu- 
dinal conductivities, previously obtained analytically for Fermi energies lying 
deep down in the lowest tail of the lowest Landau band, have been confirmed 
by our numerical investigation. However, the quadratic behaviour found at 
low frequencies for the real part of <t xx (uj) has to be contrasted with the linear 
frequency dependence that has been proposed for hopping conduction. 

Finally, a non-Drude decay of the a xx peak value with frequency, a^ x (0) — 
ct x ' x (lu) oc u>, as reported earlier for correlated disorder potentials, has been 
observed also in the presence of uncorrelated disorder potentials. A convincing 
explanation for the latter behaviour is still missing. 
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Appendix 

The iteration equations of the auxiliary quantities (required in l |18(l ) can be 
written as with Vn = H^,n+i = H N+1 N 



An+i = V N+1 R N+1 [A n + DjjRfy^DjfjRff^VN+i (19) 
Bn+i = V N+1 R^ +1 [B N - D N R N+1 D% + D%Rx +1 Dl f ]R N '+ 1 VN+i('20) 

Cn+i = V N+1 R N+1 [C N + D%R N ^ 1 D%]R N+1 V N+1 (21) 

Fn+i = V N+1 R N+1 [ f n + d n r n+i D n] r n+i^+ 1 (22) 

Gn+i = V N+1 R N+1 [G N + D]^ R^ +1 D}j]R N+1 V N+ i (23) 

Hn+i = V N+1 R%f +1 [H N + Dx R n+1 D n ) ]R n ' +1 Vn+i (24) 

d n+i = V N+1 R N+1 D N R% +1 V N+1 (25) 

D N+1 = V N+1 R N+1 D N R N+1 V N+1 (26) 

D N+1 = V N+1 R N+1 D N R N+1 V N+1 ( 27 ) 

D N+1 = V N+1 R N+1 D N R N+1 V N+1 (28) 

D N+l = V N +l R N+l D N R N+l V N+l (29) 

d n+i = V n+1 R n+1 D n R n+1 Vn+i (30) 

d n+i = V N+1 R N+1 D N R% +1 Vn+i (31) 

D N+1 = V N+1 R N+1 D N R N+1 V N+1 (32) 

-°7V+1 = ^V+l^+l^TV^Tv'+l^+l ( 33 ) 

-°iV+l = ^Ar+i^+i^TV^AT+l^V+l ( 34 ) 

D N+1 = V N+1 R N+1 D N R N+1 V N+1 ( 35 ) 

= V N+1 R+ +1 [E N + (hu>)l) R N+1 V N+1 (36) 

E N+ i = V N+1 R N+1 (E 2 N + (Rw) /) R N+1 V N+1 (37) 

E N+1 = V^ +1 R N ^ 1 \E% + (h^I^R^^VN+i (38) 

E N+1 = V N+1 R N ^ + (hu>)l) B%+ 1 V N+1 , (39) 
where the auxiliary quantities are defined by 

An = V N (Y,G N M(^) 2 Gt; + ~ 2ied ij l]x j Gj N \v N (40) 
Bn = 4 ( E [(^) 2 G„- + 2ie%I + (fiu) 2 G±] x 3 G^\ V N (41) 
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N 



N 



N 



N 



Y,{^fG + m x i G^ + x j G+ N \V N 

i,j 

N . 
i,j 

N s 
AT s 

£(m 2 g^ •'•^ ; , ; - r / ; :;'v ) r 

^(^G^G^+JViv 

AT . 

i ' 

N s 

Y,{^)G^t^Gi N \v N 

i ' 
N . 

^(MG^G^+Viv 

AT s 

^(MG^GrJ^ 

AT . 

^(^)G+ iXi Gt N )v N 

i ' 
N , 

Y^(hu)G Ni XiGih JV N 

i ' 
N s 

Yl( huJ ^ G Ni XiG iN \V N 

i ' 
N . 

^MG^G^Vjv 

N . 

^(MG^G^+Viv 

AT s 

^(MG^G^Uv 
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^ i ' 

En =vl(j2(^)G m G m \v N 

^ i ' 
^ i ' 

For the translation of the coordinates Xi — ► Xi — a one gets 



A N 




A N + Hj^- 2aD% + a 2 E 2 N + aH N 


B N 




Bn + -Hjv — -^iv + a (^Ar — ^1v) 


C N 




C N + H 2 N - 2aD% + a 2 + aH% 


Fn 




F N - aD% + 2a£>jv 5 - aL>^ - ai?^ - 


G N 




Gat + apjv + £>#)- 2a£>J? + a 2 ^w 


H N 




H n - a£>jv° + 2aD]v 4 - aD^ 1 - ai^ - 


£,1/2' 




n l/2 , rr3 


D 3/4' 






D 5 ' 

N 




D° N — aE% 


D 6 ' 

N 




D% - aEj, 


n 10/ll' 

N 




N + H N 


D 14 ' 

N 






n 15 ' 

U N 




D 1 ^ — aE N 


D 16 ' 

N 




D 1 ^ - aE% , 



with the following abbreviations 





= a(D 1 N + D% 


) 




= a(D% + D% 


) 




= aV! r (R N ' + ~ 




H N 


= aVk(RT- 




H% 


= aV} r (II%-- 


Rn)Vn 
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